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I. INTRODUCTION 

Percolation has been of tremendous interest to di- 
verse scientific disciplines for more than two decades. On 
one hand, the percolation problem is almost trivially sim- 
ple to state, on the other hand, it has an abundance of 
fascinating features. Although percolation is purely geo- 
metric in nature, it reflects many of the concepts of crit- 
ical phenomena. The percolation transition is the proto- 
type of a geometrical phase transition. In this respect, 
the role of percolation may be compared to that of the 
Ising model in conventional critical phenomena. To date, 
percolation represents one of the best studied areas of sta- 
tistical physics and yet it is a highly vivid and fascinating 
area of modern research. 

Directed percolation (DP) [|| is an anisotropic variant 
of ordinary isotropic percolation (IP) in which an effect 
or activity can percolate only in a given preferred direc- 
tion. If the preferred (longitudinal) direction is viewed as 
time, DP can be interpreted as a dynamic process. In the 
dynamic interpretation, DP represents one of the most 
prominent universality classes of nonequilibrium phase 
transitions: the DP universality class is the generic uni- 
versality class for phase transitions between an active and 
an absorbing inactive state [[| [| . 

Perhaps, DP is the simplest model resulting in branch- 
ing self-affme objects. It has many potential applications, 
including fluid flow through porous media undergravity, 
hopping conductivity in a strong electric field 0, crack 
propagation Q, and the propagation of surfaces at de- 
pinning transitions in one dimension Q. Moreover, it is 
related to epidemics with recovery || and self-organized 
critical models |9). 

In this paper we discuss multifractality |L0| in DP. The 
notion of multifractality describes the situation that an 



exhaustive characterization of the distribution of a lo- 
cal physical quantity requires the introduction of an in- 
finite set of independent critical exponents. This means 
that the distribution is not controlled by a single or sev- 
eral length scales, but rather by an infinite hierarchy of 
such. Systems in which multifractality has been observed 
include turbulence (TTJ , diffusion near fractal s |fl2| and 
polymers |ll|, electrons in disordered media poly- 
mers in disordered media [ |15|, r andom ferromagnets [p.61 , 
chaotic dissipative systems Jl7|, and heartbeat Jl8[ . 

Since the 1980s several publications have appeared 
that have addressed multifractality in IP by studying ran- 
dom resistor networks (RRNs) §0|, |l] H> |§ §§, |§ . 
A RRN is a simple bond percolation model on a d- 
dimensional hypercubic lattice where bonds between 
nearest neighboring sites are occupied with probability 
p with a resistor or empty with probability 1 — p. In a 
typical setup one has two leads positioned at two con- 
nected lattice sites x and x' at which a fixed external 
current / is injected and, respectively, extracted. It is 
well established that the distribution of bond currents in 
RRN is multifractal. This multifractality is accessible by 
experiments or simulations via the moments of t he dis - 
tribution that are [up to technical details, cf. Eq. ( 2.13 )] 
given by 



MP (x,x') 
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(1.1) 
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Here h denotes the current flowing through bond b and 
the sum is taken over all bonds. Some of these multifrac- 
tal moments are corresponding to quantities that have a 
particularly prominent role in percolation theory. For ex- 
ample, Mj°\ Mj \ Mj 2 \ and Mj 00 ^ are proportional to 
the number of bonds on the backbone (bonds which carry 
nonzero current), the total resistance, the noise (second 
cumulant of the resistance fluctuations), and the number 
of the red bonds (bonds which carry the full current). At 
the percolation transition, the moments are described by 
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a power law in the separation between the leads, 



(x,x') ~ \x 



(1.2) 



where v is the correlation length exponent of the IP uni- 
versality class. The i/j] p are referred to as multifractal ex- 
ponents. The multifractality of the current distribution 
manifests itself in the fact that the ip} F , when understood 
as a function of the index I, are not related to each other 
in a linear or affine fashion. 

To our knowledge the issue of multifractality in DP 
has not been addressed hitherto. Here, we present such 
a study. By employing our real world interpretation of 

tnman diagrams §|, ||, || ||, |t], ||, ||, ||, |3,J|, 
[m| as well as our concept of master operators p3[ , 
we find that the moments of the current distribution in 
DP are governed by a family {ifii} of multifractal ex- 
ponents analogous to {"0; IP }- For example, if the mul- 
tifractal moments are measured in the preferred time- 
like direction between two leads at x = (x± = 0, t) and 
x' = (x' ± = 0,t'), then we find that 



Mf ] (x,x') ~ (i-i') 



(1.3) 



where v\\ is the longitudinal correlation length exponent 
of the DP universality class. We calculate the family 
{ipi} to second order in an e-expansion about the upper 
critical dimension 5. 

Brief account of our results has been given in Ref . |35| . 
The aim of the present paper is to present our work in 
some detail. Its remainder is organized as follows: Sec- 
tion [n] is devoted to modeling issues. We start by briefly 
reviewing random resistor diode networks and Ohm's and 
KirchhofFs laws. Then, static noise is introduced in the 
networks and corresponding noise cumulants are defined. 
These noise cumulants are important because they are 
closely related to the multifractal moments that are not 
accessible by our methods directly. Then, we set up a 
generating function for the noise cumulants by employ- 
ing the replica trick. The replication leads to an effective 
Hamiltonian that we refine into a field theoretic Landau- 
Ginzburg- Wilson functional TL. We conclude Sec. || by 
analyzing the relevance of various coupling constants ap- 
pearing in TL. We show that the noise cumulants are 
associated with dangerously irrelevant couplings. Sec- 
tion III contains the core of our renormalization group 
analysis. Its main content is the calculation of the gen- 
crating function for the noise cumulants by employing 
renormalized field theory. It starts with determining the 
diagrammatic elements of our perturbation calculation 
upon recasting TL into the form of a dynamic functional. 
Next, we demonstrate how we incorporate the danger- 
ously irrelevant noise couplings into our diagrammatic 
calculation via operator insertions. We show that the op- 
erators associated with the noise cumulants are master 
operators. We derive the scaling behavior of the gener- 
ating functions. Then, the scaling behavior of the noise 
cumulants and multifractal moments is readily extracted 
via taking derivatives. Finally, our result for {■0/} is 



stated and several consistency checks are provided. Con- 
cluding remarks are given in Sec. [V. Technical details 
are relegated to two Appendices. 



II. THE MODEL 

A. Random resistor diode networks 

Our approach is based on a model which captures 
both, IP and DP, viz. the random resistor diode net- 
work (RDN) introduced by Redner Ql^, |§. We start 
by reviewing substantial features of the RDN. A RDN 
consists of a rf-dimensional hypercubic lattice in which 
nearest-neighbor sites are connected by a resistor, a pos- 
itive diode (conducting only in a distinguished direction), 
a negative diode (conducting only opposite to the distin- 
guished direction), or an insulator with respective prob- 
abilities p, p + , p-, and 5=1— p — p+ — P-- The three 
dimensional phase diagram (pictured as a tetrahedron 
spanned by the four probabilities) comprises a nonper- 
colating and three percolating phases. The percolating 
phases are isotropic, positively directed, or negatively di- 
rected. Between the phases there are surfaces of contin- 
uous transitions. All four phases meet along a multicrit- 
ical line, where < r := p+ = p- < 1/2 and p = p c (r). 
On the entire multicritical line, i.e., independently of r, 
one finds the scaling properties of usual isotropic perco- 
lation .{r = 0). For the crossover from IP to DP see, e.g., 
Ref. [§9i. 



Kirchhoff 's laws 



To be specific we choose n = l/yd(l, . . . , 1) as the 
preferred direction. We assume that the bonds bu -\ be- 



tween two nearest neighboring lattice sites i and j are 
directed so that b(ij) • n > 0. We suppose that the di- 
rected bonds obey the non-linear generalization of Ohm's 
law 



,v b( . . 



(2.1) 



where Vt V* is the voltage drop over the bond 

between sites j and i and denotes the current flow- 

ing from j to i. In the following we drop the subscript 
whenever there is no risk of confusion. The bond 
conductances — Sblb are equally and independently 
distributed random variables. The 7b take on the val- 
ues 1, 0{V), 9(—V), and with respective probabili- 
ties p, P-, and q. 9 denotes the Heaviside function. 
The nature of the <a will be specified below. Note that 
the diodes are idealized: under forward-bias voltage they 
behave as "ohmic" resistors whereas they are insulat- 
ing under backward-bias voltage. To prevent confusion, 
we point out that the round brackets in Eq. (2.1) con- 



tain the functional argument of the bond conductance. 
Because the bond conductance depends on the voltage 
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drop over that bond only via a Heaviside function and 
sign (Vh) = sign (4) we may write (14) = (4) . 

Assume that an external current 7 is injected at x and 
extracted at x' . It is understood that x and x' are con- 
nected. The power dissipated on the network is by defi- 
nition 



P = I[V X - V x .) 



(2.2) 



Using Ohm's law it may be expressed entirely in terms 
of the voltages as 

P = R+{x,x')- l [V x -V xl ] 2 

= Y,<nL<y£vl = p{{v}). (2.3) 



configurations. This average, however, requires a little 
caution because the resistance between sites not con- 
nected by any conducting path is infinite. Therefore, 
we will exclusively consider those sites x and x' known 
to be connected by such a path. In practice this is done 
by utilizing an indicator function X+{ x , x ') that takes the 
value one if x and x' are positively connected, i.e., if 7 
can percolate from x to x', and zero otherwise. Note 
that (x+{x,x'))c — (x-(x',x))c is nothing more than 
the usual DP correlation function. With help of the in- 
dicator function the average resistance, or more general 
the Ith moment of the resistance, can by written as 



{x+(x, x')R + (x, x') l ) c /{x+{x, x')) c 



(2.8) 



The sum is taken over all current carrying bonds (the 
backbone) between x and x' and {V} denotes the corre- 
sponding set of voltages. R+(x, x') stands for the macro- 
scopic resistance when 7 is inserted at x and withdrawn 
at x' . Similarly one defines R^(x,x') as the macro- 
scopic resistance when 7 is inserted at x' and with- 
drawn at x. The two quantities, however, are related 
by R + (x,x') = R^{x',x). From the power one obtains 
Kirchhoff's first law 



<i> 



-h 



(2.4) 



as a consequence of the variation principle 



d 



(2.5) 



The summation in Eq. (2.4) extends over the nearest 
neighbors of i and 7j is given by Ii = I (§i tX — 8i, x >)- 

Alternatively to Eq. (2.3), the power can be expressed 
in terms of the currents as 

P = R+ (as, x') I 2 =Y,Pk (4) 4 2 = P ({I}) , (2-6) 



with {7} denoting the set of currents flowing through the 



individual bonds. 



Pb - u k 
It is understood that pi, (7/, 



(T u denotes a bond resistance. 



7? = whenever ab (lb) = 



0. Kirchhoff's second law, saying that the voltage drops 
along closed loops vanish, can be stated in terms of the 
variation principle 



d 



dm 



P 



= 0, 



(2.7) 



i.e., there are no independent loop currents 7") circulat- 
ing around a complete set of independent closed loops. 

As customary in dealing with electric networks, Kirch- 
hoff's equations may be exploited to calculate the total 
resistance R + . Of course we are not primarily interested 
in any particular random configuration C of the diluted 
network but rather in an average (■■ - )c over all these 



C. Incorporation of noise 

In the usual RDN all the % are equal to one. Here, 
we are interested in a more general setup where the 
bond conductances fluctuate statically about some aver- 
age value ^. In other words we are interested in an RDN 
with static noise. This noise is modeled by distribut- 
ing the % according to some distribution function / with 
mean ^ and higher cumulants A^- 2 ^ satisfying <C <T ( . 
The condition on the cumulants is imposed to avoid un- 
physical negative conductances. The distribution func- 
tion / might for example be Gaussian. Nevertheless, our 
considerations are not limited to this particular choice. 

The noise average will be denoted by 



{•••}/=/ 



(2.9) 



and its Ith cumulant by 



(2.10) 



A=0 



Of course, both kinds of disorder, the random dilution of 
the lattice and the noise, influence the statistical prop- 
erties of the macroscopic resistance. Both are reflected 
mutually by the moments 



m iX+(x,x') {R + (x,x') 1 } 

r(0^ ^ _ \ J l_c ( 2 n ) 



M£(x,x') = 



and the cumulants 



(x+(x,x')) c 



r(l)( ^_ {x + {x,x'){R + (x lX ') l }f) c 
Cr[X ' X) - ( X+ (x,x')) c 



(2.12) 
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D. Moments of the current distribution 

Primarily we are interested in the moments of the cur- 
rent distribution defined by 



M 



f(x,x') = ( X (x,x')J2[ I j) ) /{X{x,x')) 



2l\ 



(2.13) 



(') 



There exists an intimate relationship between the M) 

and the that can be understood as a consequence 
of Cohn's theorem |40| |. We shall now briefly review this 
relationship. 

In Eq. ( |2.9| ) the noise average is defined as an average 
with respect to the distribution of the bond conductances 
Ob- Alternatively one might express the macroscopic re- 
sistance in terms of the bond resistances pb and average 
over their distribution. Of course, not only the (Jb but 
also the pb are distributed independently and equally. 
Assume that the distribution function of the deviations 
Spb = pb — ~p of the resistance of each bond from its av- 
erage p has the form 



and that 



g s (6p b ) = -h ( — 
s \ s 



lim g s (Spb) = 5 (Sp b ) 

s— *0 — — 



(2.14) 



(2.15) 



s is a variable with units of resistance which sets the 
scale of the distribution. With this form of g s , the nth 
cumulant v n of 5pb tends to zero as s™. This follows from 
the generating function c(As) of the v n : 



exp [c (As)] = {exp (X5pb)} f = j dy h (y) exp (Xsy) 

X 



exp 



(2.16) 



where v n — c n s n with c n being constants. In gen- 
eral {R + (x, x') 1 }^ depends on the entire set of cumu- 

lants {v n }. However, in the limit s — > the leading 
term is proportional to vi as we will see in a moment. 
Consider the generating function C (A) of the cumulants 

{R + (x,x') l } ( f c) , 

exp [C (A)] = Y\. d ^P b 9 s ftPk) exp x ')] ■ 

b 

(2.17) 

Expansion of the macroscopic resistance in a power series 



in the Spb leads to 

exp [C (A)] = Y[dy b h (y b ) exp 



XR° + (x,x') 



^2, s k d k R(x,x') 



(2.18) 

where RR{x,x') is the r esistance when Spb = for every 
bond b. Equation ( [2.18] ) can be rearranged as 



exp [C(X)} 



exp 



XR° + (x,x') + Xj2 E 



s k d k R+{x,x') 



k=2 b. 



Qk 



dzh ■ ■ ■ dz by 



Tl cxp c ( 
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k\ dpb ± ■ ■ ■ dpb k 



= exp 



i=i b \ ° Pb p/ 



(2.19) 



where /$ is a function of As 1 . Hence, for I > 2, 

l 



{R+^x') 1 }^^^ 



dR + (x, x') 



dpb 



1 + O(s) 



(2.20) 



In the limit s — > the leading term is 

dR+{x,x') 



{r + (x,x'Y}^ = Vl J2^- 



dpb 



b v 



2/ 



(2.21) 



where we have used Cohn's Theorem in the form 

2 



dR+(x,x') _ (\ 
dp b ~\I 



(2.22) 



Upon substitution of Eq. (2.21) into Eq. (2.12) one finds 
for the noise cumulants 



C^{x,x')=viM^>(x,x') 



(2.23) 



i.e., the ith noise cumulant is proportional to the Ith 
multifractal moment of the current distribution. In the 
following we will exploit this relation in the sense that we 
will study the as a surrogate for the Mj . We will 

see below that the are accessible in an elegant way 
by renormalized field theory whereas we do not know of 
such an approach for the Afj^ directly. 
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E. Replication 

In this section we devise a generating function for the 
C)j based on the ideas of Stephen pl| and their refine- 
ment by Park, Harris, and Lubensky (PHL) |22| . We 
demonstrate that this generating function indeed serves 
its purpose and explain how the average conductance can 
be extracted from it. 

PHL introduced the quantity 



i>j{x) = exp (iX ■ V x j , A ^ 



(2.24) 



V x is a (D x £')-fold replicated variant of the voltage V x 
at lattice site 



V, 



(i,i) 



V, 



(i,£>) 



V x = 



(2.25) 



(-5,1) 



(£,D) 



A is, apart from a factor — i, a replicated external cur- 
rent that is like V x a matrix with (D x E) compo- 
nents. The corresponding scalar product is defined as 

^■v x = Y, D J^ {a ' 0) v^ p) . 

In order to proceed towards the desired generating 
function we now consider the two-point correlation func- 
tion of ?p~(x) 



where the average is defined by 
/ f 1 



rep 



(2.26) 



x exp 




Here, dVj is an abbreviation for [] a '„ =1 dV^ 



(«5) 
6 



7=1 b 



V, 



(7,5) 



(2.28) 



is the replicated version of the electric power with 

(1,5) V (D,5) 
, .c 7 ■ ■ ■ ! vx 

in Eq. ( [2.27 ) is given by 



(V x "' i ■ ■ ■ ,V X ^'"'). The normalization factor 



Z 



/n^exp[-ip(W,{af} 



C 



(2.29) 



At this point we would like to comment on regular- 
ization issues. First, it is important to realize that the 
integrands in Eqs. ( 2.27 ) and ( 2.29 ) depend only on volt- 
age differences and hence the integrals are divergent. To 
give these integrals a well defined meaning one can intro- 
duce an additional power term ^ V 2 . Physically the 
new term corresponds to grounding each lattice site by 
a capacitor of unit capacity. The original situation can 
be restored by taking the limit of vanishing frequency, 
u; — > 0. Second, it is not guaranteed that Z stays finite 
because infinite voltage drops may occur. Thus, the limit 
limc^o Z DE is not well defined. This problem may be 
regularized by restricting the voltage variables to a finite 
interval. However, we have to bear in mind that A = 
has to be excluded properly. We take care of both points 
simultaneously by resorting to a lattice regularization. 
To be specific, we switch to voltages taking discrete val- 
ues on a (D x .E)-dimensional torus, henceforth called the 

replica space. In practice we set ■& = Adk with k being 
an (D x i?)-dimensional integer with — M < k^ a '^ < M 
and fc( Q ^> = k i - a ^mod(2M). M = d M /M is the gap 
between successive voltages and §m is the voltage cut- 
off. The continuum may be restored by taking 6m —* oo 
and A6> — > 0. By setting M = m 2 , <&m — ^o TO > an d, 
respectively, A$ = z?o/ TO , the two limits can be taken 
simultaneously via m — > oo. For the replica currents we 
set 



A = AA I , AA A9 = ir/M 



(2.30) 



where I is a (D x _E)-dimensional integer taking the same 

values as k. This choice guarantees that the complete- 
ness and orthogonality relations 

J)DE J2 CX P (** ' <0 = S X,0 mod(2MAA) ( 2 ' 31a ) 



(2M) 



and 



(2M) 



-— ^exp (*A 



6~ ~ , (2.31b) 

9 , mod(2MA0) v ' 



do hold. Equation ( 2.31 ) provides us with a Fourier 
transform in replica space. It is important to note that 
the replica space Fourier transform of ^>«(x), 



$~ (a:) = (2M)~ DE exp (*A ■ £) ^(x) 



A#0 



-DE 



(2.32) 



satisfies the condition ^« ^^(x) — 0. Hence, it can be 

identified as a Potts spin with q = (2M) DE states. 
The relation of the RDN to the Potts model that emerges 
here as a consequence of the lattice regularization is im- 
portant as well as intuitively plausible. That is because 
the purely geometric properties of the RDN are those of 
percolation, either IP or DP, depending on which sector 



6 



of the phase diagram is considered. It is a well known 
fact that the Potts model describes percolation in the 
limit q — > 1 [43)] . Note that this limit corresponds to the 
replica limit via q = (2M) DE . 

Before procee ding with the evaluation of the correla- 
tion functions ( 2.26 ) we would like to make one more 
comment on the replication procedure. The replication 
scheme employed here goes beyond the usual replica trick 
in the sense that it involves a second replication parame- 
ter, i.e., that the replicated quantities are (D x £')-tuples 
and not jus t £>-t uples. Thi s subtlety has its origin in the 
definitions (2.11) and (2.12) which require to treat the av- 
erages (• • • )c and {■ ■ - }f independently. The great ben- 
efit of the replication is to provide for the free parameter 
D which we may tune to zero J44|. In this replica limit 
the normalization denominator Z~ DE goes to one and 
hence does not depend on the distribution of the bond 
conductances anymore. Then the only remaining depen- 
dence on this distribution rests in the electric power P 
appearing in the exponential in Eq. ( 2.26| ) . In the replica 
limit, therefore, we just have to average this exponen- 
tial instead of the entire right hand side of Eq. (2.26). 
This average then provides us with an effective electric 
power or Hamiltonian which serves as vantage point for 
all further calculati ons. The effective Hamiltonian will 
be discussed in Sec. II F 



No w we come back to the role of the correlation func- 
tions ( 2.26| ) as a generating func tion. The integrations as- 
sociated with the average (2.27) are not Gaussian. How- 
ever, they can be carried out in an approximating man- 
ner by applying the saddle point method as it was done 
by Harris [}45| in the related problem of nonlinear ran- 
dom resistor networks. We extract the leading contri- 
bution to the integral stemming from the maximum of 
the integrand. This maximum is determined by the so- 
lution of Kirchhoff's equations, i.e., by the macroscopic 
resistance. The conditions under which the saddle point 
method works reliable will be outlined in the next para- 
graph. Under these conditions we find 



G 



(x,x', ^ = ( J] exp 



\/3=l 



AW 



R { f (x, x') 



fi c 
(2.33) 



up to an unimportant multiplicative constant that goes 

to one in the replica limit D ~ -> 0. On defining Ki(X) = 

Es=iEa=i('' 1 ' '' 3 ') 2 ] 1 we obtain by expanding in terms 
of cumulants 



G 



exp 



.1=1 



-1/2)' 



K x (I) (x,x') l } ( 



c 

(2.34) 



Upon expanding the exponential we get 
G (x,x', a) = (x+(x,x')) c 

x { 1 + ^klK Kl (x)C${x,xi) + - 



i=i 



(2.35) 



We learn here that C# can be calculated via 

Cf {x,x') = { X+ {x,x'))- c } 



G ( x, x , A 



A = 



(2.36) 



i.e., G represents indeed the desired generating function 
for the noise cumulants. 

Now to the conditions for the saddle point approxima- 
tion. We work near the limit when all the components 
of A are equal and continue to large imaginary values. 
Accordingly we set |45]] 



0) 



(2.37) 



with real A and , T,a=i^f > = °- Thc saddle 

point approximation gives Eq. (2.33) provided that 



I An | » 1 
On thc other hand one has 



(2.38) 



(2.39) 



Thus, one can justify the expansion in Eq. ( 2.3 5| ) by in- 
voking the conditions 



Xl^D- 1 and e 1/3)2 «l 



(2.40) 



It is important to realize that the replica limit D — > 
allows for a simultaneous fulfillment of the conditions 
( pT38[ ) and ( gig ). 



Field theoretic Hamiltonian 



This section presents the derivation of our field theo- 
retic Ham iltoni an for the noisy RDN. We start by revis- 
iting Eq. (|2.27| ) and note that the effective Hamiltonian 
heralded in Sec. [IE is given by 



In yexp 



>(«)]) 



(2.41) 



Carrying out the average over the diluted lattice config- 
urations leads to 



H,, 



(2.42) 
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where we have introduced 



sites i and j we write 



K (jpj = In < q + p Yl i os.\> 

V ' I /3=1 



2 



r r (p) i i 

+ p+n ri cx p -^—e^^)^ 2 1 

,3=1 ^ a=l L -I J 



/8=: 

E , D 

p- n n ex p 

/3=1 k Q = l 



Jfl) 



/J 

(2.43) 



Now recall our choice for A in Eq. (2.37). Since A and $ 

are related via Ohm's law Eq. ([2~l|), we have to choose $ 
consistently: 



(2.44) 



with real i? and C (q,/3) , E«=i C (q,/3) = 0- We impose the 
conditions 

|i?o|»l, Vl^D- 1 , C (/3)2 <1- (2.45) 

Under these conditions we have 6{^ a ^) = 0(& o ). 
Hence, we can write 



K\0) = lnjg+]> + 0(-i? o )+p_0(i? o ) 

+ [p + p + 0(0 o )+P-0(-0o)] 
x exp 



gHZ^A 1 "^ (7 



.1=1 



l\ 



(2.46) 



After some additional algebraic steps and by dropping a 
term 

9 (0„) In [1 - p - p+] + 6 (-i? ) In [1 - P - P-] , (2.47) 

that does not depend on the bond conductances we arrive 
at 



K 



(fj = 9 (tf ) K+ (#) + (-0 O ) if_ (tf) , (2.48) 



with 



A± (V) 



In < 1 + exp 

1 - P ~ P± 



U=l 



(2.49) 



In order to proceed towards a field theoretic Hamiltonian 
we expand K (■&) in terms of the tp^(i) = exp(iA-i?j). Re- 
calling that b is a shorthand for the bond 6 <i 3> between 



A' 



EE 



exp 



K Q 



(2M) 

E ^A (0 ^-a (i) 

A/0 

[9(X )K+(\) +9(--X,- l K (A) . (2.50) 



Here, we have exploited that 0($o) = 9(Xo). A±(A) 
stands for the Fourier transform of A±($), 



K- 



P)^(2i7r? exp 



iX • 



A± p) . (2.51) 



For evaluating Eq. ( [2.51 ) we switch back to continuous 
voltages, 



Ka 



(?) 



dd exp 



(-iX ■ Tj 



hx< 1 



x exp 



P + P± 
1 - P ~ P± 

(2.52) 



where we have dropped a factor (2§m) DE - Taylor ex- 
pansion of the logarithm yields a series of terms of the 
form 



di9 exp 



—iX ■ 7 - aid 2 -J2 b i K i (*) 

1=2 



(2.53) 



where a the 6/ are constants of order O (s°) . For nota- 
tional simplicity we dropped the subscript ±. In addition 
to the expansion of the logarithm we expand in a power 
series in s, 



-iX ■ tf-o?i?' 



Eq. (2.53) = / dtf exp 



x <M 



OO 1 

Hp) ■ (2-54) 



1=2 



Here, the Pi are homogeneous polynomials of order 21 in 
A which are a sums of terms proportional to 



n«(3)' 

i>2 



(2.55) 
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such that VJ. il 



Completing squares in the expo- 



nential in Eq. ( 2.54 ) gives 
Eq. (2.54) = exp 



pOO 




/ exp 


-aid 2 


J —oc 






(2.56) 



where we have omitted multiplicative factors decor ating 
the p. Due to the homogeneity of the p, Eq. ( 2.56 ) can 
be rearranged as 



Eq. (2.56) 



exp 



exp 



A 2 " 
4as" 


f- 


OO 

^E 

1=2 


Pl-r 


p>- 


1 


A 2 " 
4as" 


f- 


OO 

^E 

/' = ! 



-'p, 



(?) 




l + 0(s) 




up to multiplicative factors. By keeping only the leading 
contributions, we deduce for K(X) that 



E w ±-p 
P =i 



A 2p 



E r 

p, 



(2.58) 



with r±, U>± lP ~ and f±.p, ~ A^ /s 21 being expan- 
sion coefficients. ^ 

The terms w n X 2p are irrelevant in the renormaliza- 



tion group sense for p > 2, cf. Sec. |IIG| . The up,p(A) 
are irrelevant as well. However, we will demonstrate in 



Sec. [II G| that the terms proportional to Ki{\) are indis- 
pensable in studying the noise cumulants; they are dan- 
gerously irrelevant. Therefore, we restrict the expansion 
of K±(\) to 



K± (a) = t± + w±\ 2 + v ±-i R i P 



1=2 



(2.59) 



with w± = w±,i, and v±j = v± t Ki- Nevertheless, 
the neglected terms will regain some importance later 
on since they are required for the renormalization of 
the v± t i. The expansion coefficients in Eq. ( 2.59] ) 
satisfy T ± (p,p + ,p_) = r T (p,p^,p + ), w±{p,p+,p-) = 
w T (p,P-,p+), and v ±l i(p,p + ,p-) = v T ,i(p,p-,p+). The 



Ki are homogeneous polynomials of order 21. For I > 2 
they are possessing a S[0(D) E ] symmetry whereas A 2 
has a full O(DE) rotation symmetry in replica space. 

Next we decompose K('db) into two parts, one being 
even and the other being odd under A — > — A: 



A^O 



K+[X)+K-[\ 



+ 2^(Ao)-^(-A )] 



K + I A I — ( A 



(2.60) 



Then we insert Eq. (2.60) into Eq. (2.42). We also carry 
out a gradient expansion in position space which is jus- 
tified because the interaction is short ranged. We find 



A^O 



1 



K+ A + K- A 



x 

1 



i + ^fe-v) 2 



+ 2^(Ao)-^(-A )] 



i^+^A 

x ^_v(i)[& i -V + ---]VvW 



^ (0 



(2.61) 



We proceed with the usual coarse graining step and 
replace the ipx{i) °y or< i er parameter fields ^v( x ) which 

inherit the constraint A ^ 0. We model the correspond- 
ing field theoretic Hamiltonian TL in the spirit of Landau 
as a mesoscopic free energy and introduce the Landau- 
Ginzburg- Wilson type functional 



H = 



i> x (x)K V,A (x) 



A#0 



9 
6 



£ ^_ X (x)^ X ,(x)^ +x ,(x)| 



A,A',A + A'^0 



(2.62a) 



where 



X V, A = t 



V 2 + w\ 2 + v i R i (a) 



i=2 



+ [0(A o )-0(-A o )]r-V. (2.62b) 



In Eq. ( 2.62] ) we have discarded terms of higher order in 
the fields that are irrelevant in the renormalization group 
sense. The parameter r is the coarse grained relative 
of r+ + t_ . It specifies the "distance" from the critical 
surface under consideration. In mean filed theory the 
transition occurs at t = 0. w ~ is the coarse grained 
analog of w + + io_. The stem from v + j + U— The 
vector r lies in the preferred direction, r = rn. r, u>, 
and r depend on the three probabilities p, p+, and p_. 
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r is zero if p + = p_. We will see as we go along that 
our Hamiltonian Ti describes in the limits w — ► and 
vi — > the usual purely geometric DP. Indeed W leads 
for u> — > and v i — > to exactly the same perturbation 
series as obtained in fi 



47 



G. Relevance 

In this section we show that the vi are dangerously ir- 
relevant. The notion of dangerously irrelevant variables 
was coined by Fisher |Q . It applies to variables that can- 
not be taken to zero because the quantity examined cither 
vanishes or diverges in this limit. Later on the notion was 
carried over to field theory by Amit and Peliti p9[ . A 
characteristic feature of dangerously irrelevant variables 
is that corrections due to them determine the asymp- 
totic behavior of quantities with the above property, so 
that their effect is felt arbitrarily close to a phase tran- 
sition fl50| . In contrast, usual irrelevant variables cause 
corrections to scaling that vanish at criticality. 

We will see in a moment that the vi are irrelevant on 
dimensional grounds, i.e., they are associated with a neg- 
ative naive dimension. However, we cannot simply take 
the vi to zero by appealing to their irrelevance, because 
the amplitudes of the noise cumulants vanish in this limit. 

Now we embark on a scaling analysis in the current 
variable by rescaling A — > b^ 1 X. Here b denotes a scaling 
factor and should not be confused with the index labeling 
the bonds. By substituting i/'x( x ) = ^l-i^^ ' ln ^° ^ e 
Hamiltonian we get 



H 



o A 

^ A#0 

9 - V ^ «(x)V f « (x) 



b- 1 \+b 



» (x) 

-1 \, \ > 



(2.63) 



Renaming the scaled voltage variables \'=b X A leads 
to 



H 



"0« t (x);r, r, ,w,{v[} 

= J ddx u E 4t w K ( v ' w 

A^O 

A : A',A + A'^0 ' 

(2.64) 

Obviously, a scaling of the current variable results in a 
scaling of the current cutoff Xm — AAA/ = nm/'do, viz. 



Am — > b~ 1 \M- However, by taking the limit D — > and 
then 777, — > oo , the dependence of the theory on the cutoff 
drops out. In other words: Am is a redundant scaling 
variable. Thus, one can identify A^ and A which leads to 
the conclusion that 



= W[^(x);r,r,6 2 «;,{&%}] 



(2.65) 



Now we consider the consequences of Eq. ( |2.65| ) for the 
correlation functions of the field ipj>(x) defined by 



G 



n ({x, A | ;T,r,w,{u ; } 



xexp(-H[V'x(x);r,r, w,{vi}]) , (2.66) 

where Dip indicates an integration over the set of vari- 
ables {-0^(x)} for all x and A. Equation ( 2.65| ) implies 
that 

Gat fjx, A | ;T,r,w,{viYj 
= G N ({x, b- 1 1} ; r, r, &V . (2.67) 

From Eq. (2.67) in conjunction with Eq. ( 2.35| ) we deduce 
^(A)4 Z) ((x,x');r,r,u;,K.}) 



-21 



Ki 



(x)c^(M;r 7 ry W ,{b 2k v k }) 



(2.68) 



We are free to choose b 2 — w 1 . This choice gives 
Cr ((x,x');r, r,w, {v k }) 



^/i((x,x');r,r,{^}) , 



(2.69) 



where fi is a scaling function. We learn from Eq. ( 2. 69] ) 
that the coupling constants v k appear only in the com- 
bination Vk/w k . A trivial consequence of the fact that 
the Hamiltonian ( 2.62j ) must be dimensionless is that 



w\ 2 ~ /i 2 and A) ~ /i 2 , where \x is an inverse 

length scale. In other words u>A 2 and VkKk(X) have a 
naive dimension 2. Thus, Vk/w k ~ fi 2 ~ 2k and hence the 
Vk/w k have a negative naive dimension. This leads to 
the conclusion that the Wfc are irrelevant couplings. 

Though irrelevant, one must not set vi — in calculat- 
ing the noise exponents. In order to see this we expand 
the scaling function // in Eq. ( p. 69 ), 



c£ ) ((x,x');T,r 1 «;,{«*}) 



,(0 v i 



(I) Vl + 1 



l+l w l+l 



} , (2-70) 



with cjp being expansion coefficients depending on x, 
x', r, and r. It is important to recognize that = 
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because the corresponding terms are not g enerated in the 
perturbation calculation. Equation ( 2.70| ) can be rewrit- 



ten as 



C«((x,x / );r,r, «;,{«*}) 



C 



a 



(l) vt+i 



l+l 



WVl 



(2.71) 



Note that Vi + k/(w Vi) ~ p, , i.e., the corresponding- 
terms are irrelevant. The first term on the right hand 
side gives the leading behavior. Thus, vanishes upon 
setting vi = and we cannot gain any further informa- 
tion about cjp. In particular we cannot determine the 
associated noise exponent. In other words, the vi are dan- 
gerously irrelevant in investigating the critical properties 



of the G 



(;>2) 

R 



We conclude this sec tion by demonstrating that the 
w p >2 appearing in Eq. ( 2.58 ) as candidates for entering 
7i are irrelevant. Suppose that we had retained these 

terms. Each of them had contributed a term w p X 2p to the 



kernel in Eq. (2.62). From the analysis above it is evident, 
however, that w p had appeared in the noise cumulants 
only as w p /w p ~ /i 2 ~ 2/ . We conclude that keeping the 
w p >2 had produced only corrections to scaling and that 
neglecting them in studying the leading behavior at the 
critical point is indeed justified. 



III. RENORMALIZATION GROUP ANALYSIS 



In this section we calculate the generating function 
G(x, x',A) by employing field theory augmented by 
renormalization. For background on these methods we 
refer to Ref. |Q . We perform a diagrammatic perturba- 
tion calculation up to two-loop order. 



A. Diagrammatic elements 

Instead of working with Ti directly, we recast H into 
the form of a dynamic functional [^2|, [5^, |54| . This strat- 
egy is convenient because it simplifies the following cal- 
culations from the onset. Assuming r ^ we introduce 
new variables by setting 



rpt , ip — \r\ s , 



\r\ 1,2 g. 



(3.1) 



On substituting Eq. ( p\l| ) into Eq. ( 2.62| ) we obtain the 
dynamic functional 



J2v l K l Cx))+(9(X )-e(-X )) 



1=2 



at 



s» (xi,t) 




(3.2) 



where dj_ = d — 1. In Eq. (3.2) we dropped a term con- 
taining a second derivative with respect to t because it 
is irrelevant compared to the retained term containing 

a/at. 



From Eq. (3.2) we gather the diagrammatic elements 
contributing to our renormalization group improved per- 
turbation calculation. Dimensional analysis shows that 
the coupling constant ~g has the naive dimension (4 — 
d±)/2, i.e., d = d± + 1 = 5 is the upper critical dimen- 
sion. The Gaussian propagator G(x±,t, A) is determined 
by the equation of motion 



p(r- V 2 + uiA 2 ) + 



>(A o )-0(-A o ))^ 



xG (xj_,i, A I =S(x±)5(t). 



(3.3) 



Note that we have not included the terms proportional 
to the vi due to their irrelevance. The prop er treatment 
of these terms will be explained in Sec. [II B . Solving the 
equation of motion is straightforward. For the Fourier 
transform G(p, t, A) of G(x±, t, A) we obtain 



G (p, t, AJ = G+ (p, t, AJ + G_ (p, t, AJ , (3.4) 
where p is the momentum conjugate to x^ and 
G± ( P ,t, A) 

Ttp (t + p 2 + wA 2 ^ 



9 (±i) 9 (±A ) exp 



A,0 



(3.5) 



For the perturbation expansion it is sufficient to keep 
either G+(p,t, A) or G_(p, t, A) and hence we discard 
G_(p, t, A). The factor (1 — J« «) on the right hand side 

of Eq. (3.5) enforces the constraint A ^ 0. Due to this 
factor the principal propagator G+(p,t, A) decomposes 
into 

G+ (p, t, A) = G cond (p, t, A) - G ins (p, t) . (3.6) 
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FIG. 1: Decomposition of the primary two leg diagrams (bold) 
into conducting diagrams composed of conducting (light) and 
insulating (dashed) propagators. It is important to bear in 
mind that the conducting diagrams inherit their combinato- 
rial factor from their bold diagram. For example, the dia- 
grams A and B have to be calculated with the same combi- 
natorial factor i. 



The first part 

G cond (p,i,A) = 9(t)9(X ) 



x exp 



-tp(r + p 2 + w\ 2 



(3.7) 



is carrying replicated currents and hence we call it con- 
ducting. 



G ins (p,t)=d(t)exp[-tp (r + p 2 )]^ 



A.O 



(3.8) 



on the other hand is not carrying replicated currents and 
we refer to it as the insulating propagator. The decompo- 
sition of the principal propagator allows for a schematic 
decomposition of the principal diagrams into sums of con- 
ducting diagrams consisting of conducting and insulating 
propagators. In Fig. [l] we list the result of the decompo- 
sition procedure up to two-loop order. 

According to our real world interpretation ^3| the 
conducting diagrams may be viewed as being resistor 
networks themselves with conducting propagators corre- 
sponding to conductors and insulating propagators corre- 
sponding to open bonds. In our interpretation the times 
t appearing in the conducting propagators correspond to 

resistances and the replica variables A (up to a factor — i) 
to currents. Just as the physical currents are conserved 
in the nodes of real networks, the replica currents are 
conserved in each vertex and we may write for each edge 
i of a diagram, A^ = Ai(A,{«;}), where A is an external 



current and {k,} denotes a complete set of independent 
loop currents. 



B. Diagrammatic treatment of the dangerously 
irrelevant couplings 

Since the coupling constants vi are irrelevant, they can- 
not be treated in the same fashion as the other coupling 
constants r and w pertaining to the bilinear part of J . 
Suppose we would naively add —tpy\viKi(\) to the ar- 
gument of the exponential in Eq. ( |3.5|) and then use the 
resulting expression as the Gaussian propagator. Doing 
so we would ruin our perturbation expansion. This can 
be understood by expanding any of the diagrams in terms 
of vi. For increasing orders of this expansion one en- 
counters increasing orders of superficial divergence. The 
diagrammatic expansion can be fixed, however, by trun- 
cating the expansion in the V[ at linear order. This is 
equivalent to treating V[ by means of the insertion 



0(0 = _£!»..* 



x (P, (-P,*) , 



(3.9) 



that is a ssociated with the coupling constant vijw 1 . In 
Eq. (3J5), (f"(p,t) denotes the Fourier transform of 



s x (x ±) i). 

Now we analyze the structure of the conducting dia- 
grams after O"^ has been inserted into one of their con- 
ducting propagators. This situation is sketched in Fig. ||. 
Any of these diagrams has a current-dependent part of 
the form 



-Upw 



exp 



-Upw' ^2 K l (^i) ex P 



pw tj A 

3 

pwP (A, {k} 



(3.10) 



where P denotes the electric power of the diagram. The 
summation is carried out by completing the squares in 
the exponential. The corresponding shift in the loop cur- 
rents is given by the minimum of the quadratic form P 
which is determined by a variation p rinc iple completely 
analogous to the one stated in Eq. (2.7). Thus, com- 
pleting of the squares is equivalent to solving Kirchhoff 's 
equations for the diagram. It leads to 



x exp 



- pwR ({£}) A 2 - pw 2J Bi t j ({*}) Ki ■ *K 3 

(3.11) 
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A- nd = Ci({t})\ is the current induced by the external 
current into edge i. Ci({t}) and C,.j({0) are homoge- 
neous functions of the times of degree zero. Bij({t}) 
and the total resistance of the diagram R({t}) are ho- 
mogeneous functions of the times of degree one. By a 
suitable choice of the /tj the matrix constituted by the 
Bi j is rendered diagonal, i.e., B^j ~ Sij. At this stage 
it is convenient to switch to continuous currents and to 
replace the summation by an integration, 



{*} 



(3.12) 



i=l 



where L stands for the number of independent conduct- 
ing loops. This integration is Gaussian and therefore 
straightforward. In the limit D — > one obtains 



-t i( w l K x (Ij nd ) + • • • 

= -U<k {{t}f pw l Ki (I) + 



(3.13) 



The terms neglected in Eq. (3.13) are not required in 
calculating the ipi. This issue is discussed in detail in 
Sec. IIIC. Diagrammatically, the calculation scheme can 
be condensed into Fig. |[ Appendix [A] illustrates the 
calculation scheme in terms of an example. 



C. Renormalization and scaling 

Now we will consider the renormalization of the field 
and the va riou s parameters appearing in the dynamic 
functional (3.2). Most of the techniques we are going 
to use, such as dimensional regularization and minimal 
subtraction, belong to the standard repertoire of renor- 
malized field theory, cf. Ref. pH - The renormalization 
of the vi, however, is somewhat intricate. Hence, we will 
elaborate on it in this section. 

We start with those quantities that can be renormal- 
ized in a straightforward fashion. For these we use the 
same renormalizations as in Refs. p9|, |32f: 



r — ► r 
w —* w 



Z l ' 2 s , 
Z~ l Z p p, 

Zp l Z w w . 



(3.14a) 
(3.14b) 
(3.14c) 
(3.14d) 

= Z-V 2 Z-izy3 G; -Va u VV /2 , (3-14e) 



where e — 4 — d± and p is the usual inverse length scale. 
The factor G E = (4n)- d± / 2 T{l + e/2), with T denoting 
the Gamma function, is introduced for convenience. Z, 
Z Tl Z p , and Z u are the usual DP Z- factors, which can be 
found to second order in e in the literature ||, 
can be gleaned to second order in e from Refs. 

To prepare for the renormalization of the noise cou- 
plings we proceed with reviewing some general features 




of operator insertions. An operator Oi of a given naive 
dimension [Oi] inserted one time in a vertex function 
generates in new primitive divergencies corresponding to 
operators of equal or lower naive dimension. Thus, one 
needs these newly generated operators as counterterms in 
the Hamiltonian. The operators of lower naive dimension 
can be isolated by additive renormalization, 



Oi -^O l = O l 



[0 3 \<[0; 



X iyj Oj 



(3.15) 



Dimensional regularization in conjunction with minimal 
subtraction leads to X^j containing a monomial in r as 
a factor that is at least of degree one. Hence, these Xij 
vanish at the critical point. Being interested in the lead- 
ing behavior at criticality, we thus can neglect the oper- 
ators of lower naive dimension in the following. 

Now to the particular case we are interested in, viz. the 
insertion of O defined in Eq. ( |3.9[ ). Whereas the repre- 
sentation ( |3.9| ) is well suited for practical calculations we 
find it convenient to rewrite for the argumentation 
in this section as 



{l) 



1 



■ w 



P 



du;ET^(A) 



2p j - 
x <f>~ (p, u) (-p, -w) 



(3.16) 



where 0^(p, oS) stands for the Fourier transform of 
<y9«(p, t). Inserting O"' into a conducting diagram with n 
external legs, see Fig. ||, generates primitive divergencies 
which must be canceled by counter terms of the structure 



Pr A 



P 2a w Vv (P,t) r 



where 



(3.17) 



(3.18) 



with iri = r, is a homogeneous polynomial of degree 
2r. The notation that we use here and in the follow- 
ing is symbolic. Such a counter term depends in general 
on an entire set of external momenta, frequencies, and 
currents. 0«(p,o->) n is for example an abbreviation for 
YYi—i 0V (Pii^i). We drop constants, integrals etc. for 
notational simplicity. 

The naive dimension of operators of the type ( [3.17 ) is 
given at the upper critical dimension by 2(r + a + b + 
n — 3), as straightforward power counting shows. Hence, 
operators having the same naive dimension as O"^ satisfy 



l+2=r+a+b+n. 



(3.19) 



For n — 2 one is led to I > r, i.e., the insertion of 
generates operators containing homogeneous polynomials 
in the replica currents of degree equal or lower 21. In 
particular O^ generates an operator of type 



viKi 



(A) X (p,cf 



(3.20) 
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KijXi) 

Pi; A'^jt-j J,/ V- pi, Ai,ti 



-p, -A,t 




p, A,0 



Pi ) i ) J;- AT Pi s i 5 ^i 




p,o + 



FIG. 2: Calculation scheme. The hatched blobs symbolize an arbitrary number of closed conducting loops. The solid dots 
indicate insertions. 



P n , X„,u„ 




P r (A) 

p n , A n ,a> n y . . .V pj, A i , wa 



-p,-A,-w 




p, X,uj 



FIG. 3: O"' inserted into a conducting diagram with n ex- 
ternal legs. 



FIG. 4: An operator of the type in Eq. ( pM7D with n > 3 
inserted into a conducting two-leg diagram. 



The important question now is, if the other operators 
generated by gene rate O^ 1 also. Consider n > 3. 
With help of Eq. fl3.19| ) one obtains Z - 1 > r > 1. The 
second inequality is a basic feature of the summations 
over loop currents in the limit D — > 0. Bearing in mind 
that maximal homogeneous polynomials of degree I — 1 in 
A are gener ated, we reinsert these operators of the type 
in Eq. ( 3.17 ) with n > 3 into two- leg diagrams, see Fig. 
The resulting terms are of the form 

P r , (A)p 2o V<^(p,a,) 2 , (3.21) 

with the leading contributions satisfying r+a+b+n — 3 = 
r' + a' + b' - 1. Thus, r' > r + a - a' + b - b' + 1, or 
in other words, the homogeneous polynomials in A may 
have a higher degree than 21. Nevertheless, they are of 
the type 



2<i<r 



(3.22) 



with iTi < r < I — 1 and X)i* r i + s = r '- These 
polynomials have a higher symmetry than the original 
Ki. 

We have learned in the preceding paragraph that O"' 
generates itself and an entire family of new operators. 
Of course, the entire family of operators associated with 
has to be taken into account in the renormalization 
procedure, leading to a renormalization in matrix form 



6 «) 



O 



(0 



(3.23) 



The vector 



6 w = (o®,di t >,- 



(3.24) 



represents the family of operators associated with O^ 1 . 
The factor [Z w Z~ x )~ l reflects the fact that we have in- 
corporated a w l in the definition of O®. 

The important conclusion from our reasoning above 
is that the operators generated by 0® do not in turn 
generate O"'. For operators with this outstanding fea- 
ture we have introduced the notion of master opera- 
tors p^| . Master operators are associated with renor- 
malization matrices 



Z w = l + 0(u) 



(3.25) 



where 1 stands for the unit matrix, of a particularly sim- 
ple structure, 



( z® o ••• o\ 
o ••• o 



\ o ■ • ■ o / 



(3.26) 



The O symbolizes arbitrary elements. We will see in 
Sec. [II D, as a consequence of the simple structure of 
Z_"\ that the operators induced by can be neglected 
in calculating the scaling index of their master. Owing 
to their subordinate role, we refer to these operators as 
servants. 
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D. Scaling 



functions of a single flow parameter £ we write 



In this section we set up a Gell-Mann-Low renormal- 
ization group equation (RGE). Its solution will provide 
us with the scaling behavior of the order parameter cor- 
relation functions and finally with the scaling behavior 
of the C 



(0 

R ■ 



The bare (unrenormalized) theory has to be indepen- 
dent of the length scale pT 1 introduced by renormaliza- 
tion. Thus, the bare connected N point correlation func- 
tions satisfy the identity 



^G N ({x±,fc,w\ 2 y,T,f) 6m =0. (3.27) 



,(') 



The subscript O indicates th at the corresponding op- 
erator has been inserted. Eq. ( 3.27 ) translates via the 
Wilson functions defined by 



7... («) = /x^-lnZ.. 



K (it) = /i 

Cw (u) = M 

C P (u) = A* 

7 W (w) = - 



dp 
d\nr 



= (-£ + 7 + 2 7 P - 7«) u i 



dp 
d In ui 



9/i 
dlnp 



= 7p - 7r 



7p 7u 



(9// 
9// 



7-7p: 





lnZ« 



(3.28a) 

(3.28b) 
(3.28c) 

(3.28d) 

(3.28e) 

(3.28f) 



(the |o indicates that bare quantities are kept fix while 
taking the derivatives) into the RGE 



^ N 
^ + ir7 



1 + 7 W -^I 



xG N (jxj_,/9t, wA 2 | ;T,u,fj^j = 0. 



(3.29) 



Here, is a shorthand for 



d d d d d 

"o + P^- + TKj- + wCtu^ h PCp^- 

a^i aw ar aw op 



(3.30) 



To solve the RGE we employ the method of character- 
istics. Considering the ingredients of the RGE as being 



dp 
~d£ 
du 

~dl 



di 



hip 



d 
d 

I— lnu> 
at 



,1 
d£ 



InZ = 





A 1 j 


(3.31a) 


(3{u{£)) , 


u(l) = u , 


(3.31b) 


C P («(*)) ■ 


P(1) = P, 


(3.31c) 


K (u(£)) 


f(l) = T 


(3.31d) 


&» TO) , 


iD(l) = w , 


(3.31e) 


7(«(0) > 


Z(l) = l, 


(3.31f) 


-2« (u(*)) 


. 1(1)= I- 


(3.31g) 



These characteristics describe how the parameters trans- 
form if we change the momentum scale p according to 
p — > p{£) = £p. Being interested in the infrared (IR) 
behavior of the theory , we study the limit I — > 0. Ac- 
cording to Eq. ( 3.31b| ) we expect that in this IR limit 
the coupling constant u(£) flows to a stable fixed point 
u* satisfying @(u*) — 0. At this fixed point the RGE 
simplifies to 



^ + y7 



G N (^x±,pt,w A 2 1 ;r, u* , fij 



0. 



(3.32) 



where 7* is an abbreviation for j(u*), 7^* for j( l \u*), 
and so on. To proceed towards a solution of the RGE it 
is important to realize that the matrix 7W* inherits the 

simple structure of Z} 1 \ viz. 



,(0* 



/7M* O 

o 



o 



(3.33) 



\ O ••• O J 



By virtue of this structure, |1) = (1,0, ■•• ,0) T is a 
right eigenvector with eigenvalue 7W*. We denote the 
remaining right eigenvectors with eigenvalues 7^>2 by 
\k). The left eigenvectors are (1| = (1,0, ••• ,0) and 
(k\ = (0, O, • • ■ , O). Employing spectral decomposition 
we recast 7 (0* in terms of its eigenvalues and eigenvectors 
as — 



1 



(0* - 



\l)l {l) *{l\+Y.\ k )lk*(k\ . (3.34) 



fc>2 



Now we substitute the decomposition ( |3.34 ) into the 
RGE ( 3.32j ). Multiplying the resulting equation from the 
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left hand side with (II leads to 



N 



<G N (jxj_.pt, wA 2 } ;r, 



u ,fi\ = , 
'aw 



Here, A® is an abbreviation for 

fc>2 

Note that 



<fc|<2 (0 = ^0d( 





m J 



(3.35) 
(3.36) 

(3.37) 



m>2 



i.e., the RGE ( 3.35) co ntains all the information on the 
OW. In the from (|3.3E|) the RGE is readily solved. With 
help of the characteristics we obtain 

G N (jx ± ,p<, wA 2 } ;t,k,/xJ 

_ £7*AT/2+ 7 (i) *-iC 

(3.38) 

To account for t he naive dimensions of the various quan- 
tities, Eq. (3.38) needs to be supplemented by a dimen- 
sional analysis. Simple power counting reveals that 



,u,fi 



G N (\x-x,pt, wA 2 } ;r, 
xG/. N /iXj_, [i 2 pt, /i~ 2 u;A 2 j ; /j~ 2 t, u, l^j 



AC 



(3.39) 

Equation ( 3.38 ) in conjunction with Eq. ( 3. 39] ) now gives 
Gat f{x_L,pi, wA 2 | ;t,u, fij 

— g(d ± +r 1 )N/2-4 H /u ± +l<P/u ± 

xG N ({£x ± ,i z pt,e-*/ u± w\ 2 y,£- 1/u± T,u*,^^ . 

(3.40) 

Equation ( |3.40| ) features the well known critical expo- 
nents for DP that have been calculated previously to 
second order in e H p7| : 



V = 7' 



* = 2 + c; = 2 - - <! 1 



25 161 /4 
1 In - 

288 144 V 3 



2-k* 2 16 



67 59 /4 
288 + 144 n V3 



107 17 /4 

288 ^144 n V 3 



(3.41a) 

e 

(3.41b) 

■] 

(3.41c) 



= (2 — £*) is the resistance exponent for DP that 
we derived recently p3] to second order in s: 



= 1 



e 

24 



151 157 /4 
288 ~ 144 n V 3 



(3.42) 



i/ji is defined by ^ = f(2-7 (0 *)- The e expansion result 
of "0; is given below. 

From here only a few more steps are required to reveal 
the scaling behavior of the Gjj . Recall that our strat- 
egy is to derive the cjj from their generating function 

G(x, x'; A). By now, we know of the scaling behavior of 
a central ingredient to the generating function, viz. we 
know that the two-point correlation function with inser- 
tion scales at criticality as 



G- 



(jx__ - x' ± \,t-t',w*\' 



AW 



— pdx+V-'4>i/ v J-+l<l>/vx 

x G 2 (l\x ± -x' ± \,l z (t-t'), e-*/ v± w A 2 ) 



AW 



(3.43) 



where we dropped several arguments for notational sim- 
plicity. In the following we set x' ± = and t' — 0, again 
for the sake of simplicity. A further ingredient to the 
generating function is the two-point correlation function 
without insertion. Its scaling behavior can be inferred 
from a renormalization group treatment similar to that 
above. This analysis is comparatively straightforward 
(cf. Rcf. H) and gives 



G 2 (jxJ,i,u>A 2 
= f^+"G 2 U\*x\,l%lr+! u ±w\*\ 



(3.44) 



Now we put Eqs. (3.43) and (3.44) together. Recalling 
that our master operator is associated with a cou- 
pling constant vi / w l we write the generating function as 



G(|xj.|,t,A) =^+"|g 2 (V|xJ,^M-^ w A 2 ) 

oo 

w l 

1=2 

xG 2 (£\ X± \,£%i-^wX 2 ) m 



(3.45) 



We have the freedom to choose the flow parameter as 
we are pleased. The choice I = |x^| _1 and a Taylor 
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expansion of the right hand side of Eq. ( 3.45| ) lead to 



G(|xxUA) =|x ± | 1 - d -"/ / ' 



xx 



;=2 



(3.46) 



where the /s are scaling functions that vanish for van- 
ishing argument. Instead of choosing £ = |xjj -1 we can 
likewise choose £ — t~ x l z . This leads then upon Taylor 
expansion to 



xjl + wA 2 ^" 

oo 



(3.47) 

with the hs being scaling functions that tend to constants 
for vanishing arguments. v\\ is defined by v\\ = v±_z. 

Possessing of the generating function we solely need 
to take the appropriate derivatives to extract the scaling 
behavior of the that is by virtue of Eq. (2.23), up 



to unimportant co nstan ts, identical to that of the Mj 
With help of Eq. @ 36| ) we deduce that 



M 



(0 



(3.48) 



if measured along the preferred direction. For measure- 
ments in other directions it is appropriate to choose a 
length scale L and to express the longitudinal and the 
transverse coordinates in terms of L: |xj_| ~ L and 
x\\ ~ L z ~ T. With this choice the scaling functions 
/ reduce to constants and we obtain 



(3.49) 



We still owe our result for the multifractal exponents. 
Since we only need to compute a single element of each of 
the renormalization matrices Z} 1 \ viz. Z^ l \ we manage 
to calculate the ipi to two-loop order. In e-expansion, our 
result reads 



1 



3 • 2 2i +! 



a(0 



6(0 In ( - 



0(e 6 



(3.50) 



The a(l) and b(l) are /-dependent coefficients taking on 
the values listed in Table Q. tpQ an d i^i s tem f rom ex- 
tending the sum over I in the Hamiltonian (|2.62 ) so that 
it comprises I = and 1 = 1. Figure |5| depicts the 



2" 

1.8 
^1-6 a 
1.4 
1.2* 

1 ^ 







2 3 

/ 



FIG. 5: Dependence of ipi on I for d± = 3 (stars), = 2 
(triangles), and dj_ = 1 (squares). 



dependence of ipi on I for e = 1,2,3 corresponding to 
d± = 3,2,l. 

We point out that our result ( 3.50| ) fulfills several con- 
sistency checks, ipo is related to the fractal dimension 
D B of DP clust ers via D B = l+ip /v± - z (cf. ^ |l|). 
Equation (3.50) is in agreement with the e-expansions 
of v u 2 % 0, and L> B (see Refs. 0) to second 
order in e. ^ is in conformity with our result for the 
resistance exponent 4> given in Refs. [^9[ [52|. This has 

to be the case because = M^. Multifractal expo- 



nents like the ipi have the general feature that they are 
convex monotonically decreasing if being understood as 
a function of the index I [|55| . Our result possesses of this 
feature. Moreover, it tends to unity for large I as one ex- 
pects from the relation of ip^ to the fractal dimension of 
the singly connected (red) bonds, d rc d = 1 + V'oo/^x — z, 
see Refs. [gij |lj 



IV. CONCLUDING REMARKS 

In summary, we derived a field theoretic Hamiltonian 
for RDN that captures the multifractality of the current 
distribution in these networks. To characterize the cur- 
rent distribution, we determined the scaling behavior of 
its moments. Each moment is governed by an indepen- 
dent critical exponent, i.e., these exponents are not re- 
lated to each other in a linear or affine fashion, as com- 
monly occurs in critical phenomena under the name of 
gap scaling. We determined the family of multifractal 
moments to two-loop order. 

Our approach thrived on two cornerstones, viz. our 
real-world interpretation of Feynman diagrams and our 
concept of master operators. The real-world interpreta- 
tion remedies the apparent complexity of the field the- 
ory. It makes the field theory more intuitive and pro- 
vides practical guidance for the diagrammatic calcula- 
tions. Being interested in some quantity in real networks 
one basically just has to determine its counterpart in the 
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TABLE I: The coefficients a(l) and b(l) appearing in Eq. 3.50 







Ibl 



2 

B83B7 



3551703259 



(i 



> 7 



o(i) 



sea 



£12. 



7077SSS 



24a 



13730Jj) 



ii2i 



0.005 > a > 



6(0 



i i 324fi20S 



0.015 > b > 



Feynman diagrams. In the present case we determined 
the multifractal moments of the diagrams to study the 
multifractal moments in physical RDN. Without the con- 
cept of master operators the renormalization group anal- 
ysis presented in this paper is hardly feasible. Since the 
multifractal moments correspond in the field theoretic 
formulation to dangerously irrelevant operators O"' they 
generate under renormalization a myriad of other irrel- 
evant operators. All these must be taken into account 
in the renormalization group. Thus, on has in princi- 
ple to compute and diagonalize renormalization matrices 
that are giants for large /. Already handling the full 
renormalization matrix associated with I = 2 to one-loop 
order is tedious. The effort is comparable to that of de- 
termining corrections to scaling associated with a (A 2 ) 2 
term in the field theory of RRN |57j] . Due to the master 
operator property of the O"' , however, it is sufficient for 
our purposes to calculate a single element of the renor- 
malization matrix pertaining to each O^ 1 , and we can 
work to two-loop order with reasonable effort. To date, 
the concept of master operators has proved to be pow- 
erful in studying multifractality in RRN and RDN. We 
expect, though, that it has many more applications. It 
might be the case, that any multifractal quantity can be 
associated in the field theoretic framework with master 
operators. This is a speculation, but it is not implausible 
at all. For example, preliminary studies of the random 
field Ising model indicate the applicability of the master 
operator concept. 

To our knowledge, the issue of multifractality in DP 
has not been studied hitherto. In particular, we do not 
know of any other work, theoretical, simulational or ex- 
perimental, that provides results suitable for comparison 
to ours. For the RRN, in contrast, the multifractal expo- 
nents of the current distribution have been determined 
by Monte Carlo simulations 0, [58) . It would be very in- 
teresting to have corresponding numerical estimates for 
the RDN. According work is in progress [^9| and will be 
reported in the near future. 



APPENDIX A: CALCULATION SCHEME FOR 
DIAGRAMS WITH INSERTION 



In this Appen dix w e illustrate our calculation scheme 
outlined in Sec. 1IIB in terms of an example. For the 



sake of simplicity we consider the simplest conducting 
diagram comprising a closed loop of conducting propaga- 
tors, namely diagram A introduced in Fig. [|. With 
inserted successively in both conducting propagators the 
mathematical expression for that diagram reads 



A 



0(2) 



x expj-pt It + k 2 + (k - p) 2 + w*k 2 + w(k - A) 2 i 
x {-)pwH\Ki (k) + Ki (X - } , (Al) 



where J k is an abbreviation for (2ir) d± J d d± k. Note 

that P(A,«T) = — t*K, 2 — t(*n — A) 2 corresponds to the 
electric power of diagram A. For practical purposes we 
switch to continuous loop currents. This step is justi- 
fied at this stage because the constraint A ^ is safely 
implemented via the decomposition of the corresponding 
bold diagram into its conducting diagrams. We obtain 



A e ,(2) = / dttexp ( —i—t 



'< J exp 



2r + k 2 + (k - p) 2 



wP ( A, K 



cCk exp 



x w < Ki (k) + Ki I X — k) }, 



(A2) 



Acknowledgments 

This work has been supported by the Sonder- 
forschungsbereich 237 "Unordnung und grofie Fluktua- 
tionen" of the Deutsche Forschungsgemeinschaft. Also, 
O. S. acknowledges support by the Emmy-Noether- 
Programm of the Deutsche Forschungsgemeinschaft. 



where we have also modified the integration variable t. 
The integration over the loop current can be simplified by 
completing the squares in the exponential. We look for 
the minimum of the quadratic form P(A, k). The min- 
imum is determined by a variation principle completely 
analogous to the one stated in Eq. ( p.7| ). Thus, complet- 
ing the squares is equivalent to solving Kirchhoff 's equa- 
tions for the diagram. After carrying out the straightfor- 
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ward momentum integration we then find 



pa 



2 (4tt)^/ 2 J 

, ILJ „ 1 . 

-*| — + 2t+-p' 
P 2 



dtt{2t)- d ±l 2 



A (2) : 

x exp 
x exp 



—R(t)wX 2 / dtn exp [— twK 21 \ 



(A3) 



The example considered here highlights two points. First, 
not only primitive divergencies proportional to K 2 {X), 

but also proportional to tA 2 , wA 2 , p 2 A 2 and (A 2 ) 2 are 
generated. Second, the basic task in computing propor- 
tional to Ki(\) is to determine the currents induced by 
the external current in the conducting propagators. 



APPENDIX B: 



DETAILS ON THE TWO-LOOP 
DIAGRAMS 



R(t) = t/2 is the total resistance of diagram A. A/2, 

and, respectively — A/2, are the currents induced in the 
conducting propagators by the external current. Upon 
integrating out the loop current we get in the replica 
limit D -> 



An(2) — 



pg 



l 



8 (4tt)^/ 2 J 



dtt 1 -^ 2 




^(aI + ^a 2 



(A4) 



Next we expand the exponential function. Then we carry 
out the remaining integration. Upon discarding conver- 
gent terms that are not required for renormalization pur- 
poses we obtain in e-expansion 



A 



e>(2) 



4— + 8r + 2(p 2 +wA 2 
P 



(A5) 



Here we sketch the computation of conducting two- 
loop diagrams with insertions. We restrict ourselves to a 
few examples. The techniques presented for these exam- 
ples can then straightforwardly be adapted to the remain- 
ing diagrams. For briefness, we will exclusively consider 
those parts of the diagrams proportional to K / ( A ) . More- 
over, we set external momenta and frequencies equal to 
zero. 

At first we consider diagram H. We start by determin- 
ing the currents flowi ng th rough the conducting propaga- 
tors. Kirchhoff's law (2^) applies to the 4 vertices of the 
diagram. This allows us to eliminate 3 of the 5 unknown 
currents (one of the vertices is inactive with respect to 
this purpose since the external current A must be con- 
served). The potential drop around closed loops is zero. 
Hence we can eliminate the two rem aining unknown cur- 
rents via the variation principle (2.7) and express all cur- 
rents flowing through conducting propagators in terms 

of the times and A. The momentum integrations are 
straightforward. They can be done by using the saddle 
point method which works exact here since the momen- 
tum dependence is purely quadratic. After carrying out 
the momentum integration we have 



H 



-w l K, 



P9 4 1 
2 (4tt)^ 



dhdt 2 dt 3 exp [-t (2ti + 2< 2 + 3* 3 )] w it. 



(2*3)^/ 2 [2t 1 + 2t 2 + |t 3 l 1+ ^ /2 i ( ' 1+<2) 



h + h + h 

2*i + 2i 2 + |*3 



+ (*l+*2+* 3 ) 



*1 + *2 + 



2*1 + 2* 2 + |t 3 



2*3 



1 *i + * 2 + * 3 

2 2ti + 2*2 + |*3 



(Bl) 



r 



where n — 21. Upon doing a little algebra we rewrite where we have introduced the abbreviations 
Eq. @ as 



H 



e><<) 



-w l K, 



27i + 2I 2 + 2 1 - n h + h 



}• 

(B2) 



1 f°° 

h = 7T ^ T - / dtxdt 2 dH exp [-t (2*! + 2* 2 + 3t 3 )] 
h [*i + * 2 + t a ] n 



(2* 3 )^/ 2 [2*i + 2*2 + |* 3 ] 



3 in+dx/2 ' 



(B3) 
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i r°° 

J-—T / dhdt 2 dt 3 exp [-r (2ti + 2t 2 + 3t 3 )] 
(«) 1 Jo 

tl^l+^ + its]" (B4) 



3 1 n+d_ L /2 ' 



(2i 3 )<W 2 [2ti + 2t 2 + §t 3 ] 

i r 00 

J 3 = TT-w- / dhdt 2 dt 3 exp [-r (2ti + 2t 2 + 3t 3 ) 

*3 [*i + h + h 

(2i 3 )<W 2 [2ti + 2t 2 + §t 3 ] 



3 -in+dj./2 ' 



(B5) 



then gives the result 



3-2-™ G\ _ A 4 2n o _ . . 

7l = ^2lTT T -7 + 37- 3 + F - M 



2F 1;3 (n) 



n+1 n+1 \3 



n+1 



Here, we have used the shorthand 



(B9) 



1 f°° 

J 4 = TTu- / dt 1 dt 2 dt 3 exp [-r (2£i + 2t 2 + 3t 3 )] 



(2i 3 )^/ 2 [2ti + 2t 2 + f t 3 ] 



3 1 n+d_ L /2 ' 



(B6) 



Now consider 7i. The integrations can be simplified by 
changing variables: t\ — > \ty, i 2 — > \t(l — x — y) and 
i 3 — ► |tx. This gives after doing the integration over i 
and y 



/i = 



2-" ^3\^ /2 r(4-d ± ) 
/ (4tt)<^ 



24 V 4 



Jo 



r^ 4 / dxx^/ 2 



x (l-x) 2 (l+x) d ^\l + -xj . (B7) 

The remaining integral over x may be simplified by sepa- 
rating its divergent and convergent contributions via Tay- 
lor expansion: 



2-" /3\^ r(4-^) rf ,_ 4 

1 " 24 \a) (4^ J 

+ x 1 - d -/ 2 (-2 + d ± -4+|) 



-d ± /2 



<■/■ (\ • i.r 



1 + 2a; - -x 



(B8) 



F ro ,i(n) = 



k—r, 



k J k — m + 1 



(BIO) 



J 2 , J 3 , and ^4 can be evaluated in the same fashion as I\. 
Thus, we merely state the results: 



3 • 2-" G 2 



4 2n 



3 



2Fi,_ 3 (n) + 



_j3_ /4\" +1 

n+1 n + 1 1 3 ) 



— '! -H! 



16 e e 



+ 



n+1 n+ 1 V 3 



n+1 



+ ln 



9— n r~<2 I 9 

/ 4 = — r ^T-A-+F lt _ 3 (n) 
16 e e 



+ ■ 



n+1 



n+1' n+1 V3/ +^UJ^ (B13) 



(Bll) 



(B12) 



Carrying out the integration and expansion for small e Upon collecting, we obtain for diagram H the final result 



J 



H, 



om = -w 



pg 4 3 ■ 2-" G\ 
2 128 e 
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^--6 + F 2i3 (n)+F 2: _ 3 (n) 



- 2Fl ^-^'- 3 + ^TTl3j + n+l V 3 



n+1 



1 



n+1 



- - In 



n+1 3 



+ 2 1 



T7- + «*i,s(n) + — 7 - — 77 , 
3e 3 n+1 n + 1 \3 



n+1 



+ -z In U 



(B14) 



r 



As a further example we now treat diagram C. We em- ploy once more our calculation scheme and determine for 
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each conducting propagator the induced external current. 
This provides us with the noise cumulants of the diagram 
and leads to 



C om = -w l Ki (AJ pg 4 {l 5 + 2J 6 + 2I 7 } , (B15) 
where we have used the abbreviations 

J 5 = TjAd- / dt idt 2 dt 3 exp [-r (2tx + 2< 2 + 3t 3 )] 

x Dgx + *2 + t 3 ]" (B16) 

[4(t 1+ t 3 )(i2+t3)-tI] rl+ ^ /2 ' 

1 f 00 

J 6 = 71 — vT~ / dtxdt 2 dt 3 exp [-r (2ti + 2i 2 + 3t a ) 
*i [2tit 2 + 2t 1 t 3 + 3t 2 t 3 + 2^] 2 

x [4(t 1+ t 3 )(t 2 +i 3 )-tir +<ii/2 ' ( ' > " 

i r 00 

J 7 = u \d \ dtidt 2 dt 3 exp[-T(2ti + 2t 2 
(fi + fa) [gtitg + Mitg + f 2 * 3 + *§] " 

[4(^+t 3 )fe+t 3 )-ir i+d " /2 



3ta 



(B18) 



Exemplarily we drill into the integral I§. It can be sim- 
plified by changing variables according to t\ — > t(x — 1), 
^2 — > i(y — 1), and i 3 — > t. Upon carrying out the inte- 
gration over t we obtain 



r(4-d ± ) . d ,- 4 

(47r) d -L 



dxdy 



l 

d t -4 



[2a; + 2y-l] ai " 4 [x + y-l] r 



[4xy - 1] 



n+dj_/ 2 



(B19) 



A further simplification can be achieved by rearranging 
the remaining integrations as 



r(4-d ± 

[x + y- 



dx 



dy 



[x + y-l] 



\xy - 



4 J 



(B20) 



We learn from Eq. (B2C) that is convergent for d± < 4. 
Hence it is legitimate to evaluate it directly at d± = 4. In 
contrast to the integrals constituting diagram H we were 
not able to evaluate I§ for arbitrary n. The technical 
difficulty is the binomial appearing in the denominator 
of the integrand of I5. For n not too large, however, the 
number of terms of this binomial (2" +2 at d± = 4) is 
manageable and one can at least carry out the integra- 
tions for each reasonable n separately. We refrain from 
stating all the results because this would be rather space 
consuming. We annotate that I§ and J7 can be treated 
in a similar fashion as ^5, except that these calculations 
are somewhat more tedious. Iq and J7 are not convergent 
like is so that in practice one has to separate divergent 
and convergent contributions as it was demonstrated in 
considering diagram H. 
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